Decay rjh — > J/ipJ/ip in light cone formalism 



V.V. Braguta^B and Kartvclishvili, V. 2 -Q 

1 Institute for High Energy Physics, Protvino, Russia 
'^Lancaster University, Lancaster, UK 

The decays of pseudo-scalar bottomonium rji, into a pair of vector charmonia, J/ipJ/ij), J/iptp', ip'tp' 
are considered in the light cone formalism. Relativistic and leading logarithmic radiative correc- 
tions to the amplitudes of these processes are resummed. It is shown that the small value for the 
branching ratio of the decay rji, — > J/ipJ/ip obtained within the leading order nonrelativistic QCD 
is a consequence of a fine-tuning between certain parameters, which is broken when relativistic and 
leading logarithmic radiative corrections are taken into account. As a result, the branching ratio 
obtained in this paper is enhanced by an order of magnitude. 

Q\ ' PACS numbers: 12.38.-t, 12.38.Bx, 

O 

O . I. INTRODUCTION 

Ever since the discovery of the T meson, there have been numerous attempts of observing the lightest pseudoscalar 
bottomonium state, 7^ . However, only recently the first experimental evidence of the existence of this meson was found 
\Q ' by BaBar collaboration, in the radiative decay Y(3S) — ■> r)b + 7 0. Its mass was found to be m rib = 9388l2';!j(,stai) ± 
2.7(syst) MeV, but our knowledge of its other properties remains rather poor. 

In [2| it was proposed to look for the rjb meson in the decay rjb — > J/i/jJ/i/j, but, despite its clean signature, this 
' pH '[ process may be hard to observe due to its extremely small branching ratio: contrary to other similar processes, 
such as the decays \b — * J/ipJ/ip 0, the rate of the decay rj b — > J/ipJ/ip vanishes at the leading order of both 
relative velocity and \/M Vb expansions. The calculations made within nonrelativisitic QCD (NRQCD) 0] yield 
Br(?]b — > J/ipJ/ip) ~ 1CP 8 — [1, Q, however in 0] it was shown that the account of final-state interaction effects 
(— 1 ' can enhance it up to about 1CP 5 . 

—i, A similar conclusion can be drawn from the comparison of the decays r/b —> J/ipJ/ip, J/ipip' , tp'ip' and the processes 
of double charmonia production at B-factories. It is now clear that these processes are greatly effected by radiative and 
relativistic corrections @, S, EE [ill, El EE 0, E EE E3, EE EE H3| • With the mass of r\ h being so close to the energy at 
which B-factories operate, it is natural to expect that the same is true for the decays r/b — > J/V'J/V', J/iprf'i ^P'^P'i an d 

(■ — I hence the consideration of these processes without accounting for radiative and relativistic corrections is unreliable. 

£ — , . This was also confirmed by the calculation of radiative corrections within NRQCD, performed in [6j|. 

CN| ' In this paper, the processes r\b — * J/^'J/^i J / ipip' , ip' ip' are considered within the light cone (LC) formalism [2lT |. 

^ I , In this approach, the amplitudes of these processes are expanded in (M ce /M b i) 2 ~ 0.1, which is sufficiently small for 

■ the applicability of the method [22| . 

Q\ \ In the LC formalism, the amplitude of a process under study is decomposed into the perturbative part, dealing 

■ with the production of quarks and gluons at small distances, and the large-distance part describing the hadronization 
of the partons. For hard exclusive processes, the latter can be parameterized by the process-independent distribution 
amplitudes (DA), which can be considered as hadrons' wave functions at lightlikc separations between the partons 
inside the hadron. It should be noted that DAs contain information about the structure of mesons and effectively 
resum relativistic corrections to the amplitude. Moreover, using renormalization group evolution of DAs, one can take 
into account the leading logarithmic radiative corrections to the amplitude. 

This paper is organized as follows. In the next section DAs for charmonium are defined, and various models for 
these DAs are discussed. In the third section, the amplitude of the decay of rjb into two vector mesons is derived. 
Finally, in the last section the numerical results and their uncertainties arc presented and discussed. 
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II. DISTRIBUTION AMPLITUDES FOR CHARMONIUM 



The amplitude of the process i% — * V1V2, with V± t 2 standing for either J/ip or ip' , can be parameterized with a 
single formfactor F: 

M = Fe^ a y iP v 2 ele P 2, (1) 
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where p±,p2 and £1,62 are the momenta and polarization vectors of V\ and V2 respectively. Hence, the width of the 
decay % — > V1V2 can be written in the form 



4tt 



(2) 

V2 , the width 



where p is the 3-momentum of a final meson in the r/t rest frame. If the final mesons are identical, V\ 
T should be divided by 2!. 

In the LC formalism, the amplitude of a hard exclusive process is expanded in the inverse powers of the hard energy 
scale Eh , which for the decay 775 — > V\ V% can be identified as M Vb . The leading order contribution in this expansion 
requires the two vector mesons to be produced with polarizations Ai = A2 = 21|, but in this case the aplitudc 
([I]) vanishes. In order to obtain a non-zero result, both vector mesons need to be transversely polarized, which in 
turn means that the helicities of the quarks in both mesons must be flipped twice, and hence leads to a suppression 
factor ~ l/(M Vb ) 2 Therefore, the decay r/b — > V1V2 is a next-to-next-to-leading (NNLO) twist process, and in 
order for the calculations to be consistent one needs DAs up to twist-4. In general, twist-4 DAs should contain terms 
corresponding to higher Fock states in addition to the "valence" charm quark-antiquark state, but we expect such 
higher states in charmonium to be suppressed, and in the following we will neglect their contribution. 

The DAs for a vector meson V with momentum p and polarization vector e can be defined as follows [23j : 



(V(p,e)\c(xh P [x,-x]c(~x)\0) = f v M v 



(ex) 



Pp 



(px) 
(px) 



1 _(ez) 



Ml 



(V(p,e)\c(x)a p x[x,-x]c(-x)\0) = f T (fi) 

+ (p p x\ - p\x p 



(px) 

(e p p\ - e\p p 
(ex) 



1 



M 2 T 2 



(px) 



2 K 



If _ \Ml 

2 [€pXX exXp > (px) 
fA(n)e p \ a pe x p a x 13 
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M 2 t 2 



(3) 



ifs{p){tx) I d^ x ^ 2 (i^) 



(V(p, e)\c(x)~fpj 5 [x, -x}c(-x)\0) 
(V(p,e)\c(x)[x,-x}c(-x)\0) 



Here [x, —x] is the gluon string which makes the matrix clement gauge invariant, £ is a dimensionless variable describing 
the relative motion of the charmed quark and antiquark inside the meson, \i is the energy scale at which the DAs are 
defined, while the constants fy and /t(a*) are defined by 



(V(p,e)|c(0)7„c(0)|0) = fyMye^ 
(V^(p,e)|c(0K„ C (0)|0) = f T (v)(e^-e uP p). 

The constants /a(a*), /s(a*) can be expressed through fy,fr as follows: 



(4) 



Mm) = 



Mm) 



1 



fv - h(n) 
/t(m) - fv 



2m c (n) 



My 

2m c (fi) \ 
My ' 



My, 



V' 



(5) 



where m c (/x) is the running mass of the c quark. 

Eqs. ([3]) contain 10 independent DAs, but only 4 of these are relevant for the calculation of the 77^ — ► V1V2 decay 
rate: (£), Y i (£), $1 ( £) and $2(£) (see below). For the first two, </?i(0 and Xi(£), we will use models proposed in 
[24l [25l 12a. |27|. In [20| it was shown that, if the higher Fock states are ignored, the functions <I?i (^) and (pa(£) can be 
unambiguously determined from the equations of motion. The same is true for the functions $2(£) and X3(£). 
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In the remainder of this section, a relation between 3>2(£)>X3(0 an( i <P \(£) i Xi (£) will be derived. The functions 
$2(0 an d X3(0 can be expanded into a series of Gegenbauer polynomials [231 ] : 



X3(a;,A*) 



E 

n=2,4. 



c„( M )Cy 2 (2a;-l) 



$ 2 (x, M ) = -(l-£ 



2 \ 



1+ E rfn(^)C I 3 / 2 (2^-l) 
n=2,4.. 



(6) 



The coefficients c„(/i) and d n (fi) are related to the moments of the functions ¥>i(£), Xi(0 through the equations of 
motion 1231. 



■<r> x = r)T + ^V^(i-^))(r- 2 )*, 



2 o /x « ■ 2 

(n+i)(i-5(/*))(r>» = (D x -%)(r)L, 



(7) 



where (£™)l,t,x,* denote the moments of the DAs </>i(£), Xi(£)> X3(£)i ^2(£) respectively, while ^f^) = 
2fv / fT(fj)(fn c (fj,)/Mv)- By solving eqs. ((7|) recursively, one can determine the functions $2 (?) andx3(£)- l n H it was 
shown, that there is a fine-tuning of the coefficients of the Gegenbauer expansion at the scale \i ~ m c = m c (/i = m c ). 
Without this fine-tuning the DAs of a nonrelativistic system would show an unphysical relativistic tail already at 
the scale /j, ~ m c . In order to get rid of this tail in the DAs $2(£) and Xa(£)> fine-tuning is required between the 
coefficients c n ,d n and the parameter 6, which is related to the wave functions <^i(£),Xi(£) [2l| : 



S(m c 



S-l (i-ga)a yi(^»M~^c) 



(8) 



III. THE AMPLITUDE OF THE PROCESS 77b -> V1V2 

The diagrams that contribute to the amplitude of the process under study at the leading order in the a s expansion 
are shown in Fig. 1. The procedure of calculating the amplitude is described in detail in 21]. This is a lengthy but 





FIG. 1: The diagrams contributing to the amplitude of the process 776 — * J/ipJ/ip at the leading order in a s 
straightforward exercise, yielding a result which looks remarkably simple: 
F = J d&d^^a^f/vi/^M^ 

+ /si(/-0/t2(a*)xi(£2,m)$2(£i,m) + /s2(m)/ti(m)Xi(£i!M)$2(&,/j) 
Here the function i?(£i, £2, /i) represents the hard part of the amplitude, 

#(£i,£ 2 ,m) = ^ fvt~ 



27 



M| b a-^)(i -e 2 )(i+ 66)' 



with the decay constant defined by 



(0|6(0) 7p7 56(0)|r ?6 (p)) = i/^p p . 



(9) 



(10) 



(11) 



At this point, some comments are in order. 
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In eq. ([9]) there is a clear separation of large- and small-distance contributions. While if (£1, £2, /•*) describes 
the hard part of the amplitude, the large-distance part is parameterized by the combination of the DAs, which 
effectively include resummation of the relativistic corrections to the amplitude. A discussion of this point can 
be found in SHI. 



2. In eq. © the dependence of the hard part of the amplitude, the constants and the DAs on the scale n is explicitly 
shown. If the process in question were a leading-twist process, one could perform an exact resummmation of 
all leading-twist radiative corrections to the amplitude, ~ a s log(M 2 b /Mj,^), simply by taking fi ~ M nh pll ]. 
Indeed, for a leading-twist process, one would use the axial gauge, in which double-logarithmic and logarithmic 
corrections only appear in the self-energy diagrams and re-scattering of final particles. The double- logarithmic 
corrections are cancelled since final particles are colorless objects, while the logarithmic corrections lead to 
the renormalization of the DAs themselves. Although the decay iji, — > V1V2 is a next-to-next-to- leading- twist 
process, all the arguments given above still seem to be applicable. Note also that in eq. © there is no divergence 
in the end-point region, |£| ~ 1, indicating that all logarithms are collected. These arguments allow us to believe 
that eq. ^ includes the exact resummation of leading logarithmic radiative corrections to all loops. 

3. Whenever NRQCD and LC approaches are used to describe the same process, one should expect some kind 
of duality between the two results. For the process — > VV this duality can be checked at the leading-order 
approximation in relative velocity of the c-quark-antiquark pair inside charmonia. In particular, by taking 
infinitely narrow DAs and the constants fx, fv and masses My, 2m c at the next-to-leading order approximation 
in relative velocity [28j . 



fv 3 ' 

My <t^> 
2m r 2 ' 



(12) 



and by neglecting all radiative corrections, one gets from eq. © 

256ttV 1 



F 



81 mf 



(13) 



which coincides with the result obtained in [5j. In these formulae, (v 2 ) is the NRQCD matrix element, defined 



as 



| X +(q £ )(D)V|y( e )) 
(0\ X +(a^\V(e)) ■ 



(14) 



As noted in [20, 122], the duality between NRQCD and LC allows us to estimate the size of power corrections. 
The idea is that if one expands the NRQCD result in powers of l/M Vb , than the first term coincides with the LC 
prediction and the second term gives an estimate of power corrections to the LC result. Thus, power corrections 
to the amplitude of the rji, — > VV decay can be estimated as ~ 4v My/Mz. 



Now we have all the ingredients needed to calculate the rates of the decays % — » V1V2. 



IV. NUMERICAL RESULTS AND DISCUSSION 



A. Input parameters 



In order to obtain numerical results for the branching ratios of the decays rjb — > J/ipJ/ip, J/ipip', tp'ip' the following 
input parameters were used: 



1. The strong coupling constant a s (p) is taken at the one loop, 

4tt 



A,log( M 7A 2 )' 



(15) 



with A = 0.2 GeV, /3 = 25/3. 
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2. The mass of the c-quark in A/5 scheme, m c = 1.2 GeV. 

3. The leptonic decay constants of the J/ip and ip' mesons fy^, fy were determined directly from experimental 
data, while the constants and fj, were calculated within NRQCD in [28| : 

{fv^f = °- 173 ± °- 004 Gev2 > (/y) 2 = °- 092 ± °- 002 Gev2 > 

(f^(Mj^)) 2 = 0.144 ±0.016 GeV 2 , (//(Af. //V ,)) 2 = 0.068 ± 0.022 GeV 2 . (16) 

4. We assume that the total decay width of the 77^ meson ^toti^b) can be approximated by its two-gluon decay 
width r(ryb — > gg) which, at the leading order in relative velocity and a s , is equal to 

r t otM = nm - 99) = Y^-ji • ( 17 ) 

5. The leading twist DAs needed for the calculations are taken from models developed in [24l. I25I [26l. [27| . 

B. Estimation of uncertainties 

The most important uncertainties come from the following sources: 

1. Model- dependence of the DAs. These uncertainties can be estimated by varying the parameters of these models 
(see [3, [2^, [2^, [23] for more details). The calculations show that for the processes r)b — > J/ipJ/ip, J/ijjip 1 , ip'ijj' 
these uncertainties are no larger than ~ 5%, 13%, 30%, respectively. In fact, these uncertainties are expected 
to be rather low, due to the property that the precision of any DA model improves with evolution [24| . 

2. Radiative corrections. Within the approach used in this paper, the leading logarithmic radiative corrections 
due to the evolution of the DAs and the strong coupling constant were effectively resummed. Although 
we argued above that this is also true for all leading logarithmic radiative corrections, there is no strict 
proof of this statement. For this reason, we estimate the uncertainty due to the radiative corrections as 
~ a s (M nb /2) log(M 2 b /(4M 2 /v ,)) ~ 50%. 

3. Power corrections. As mentioned above, this source of uncertainty can be estimated as ~ A(v 2 )M v /M^ b , which 
is the largest for the decay r/ b — > ip'tp', reaching ~ 4(w 2 )^-M|,/M 2 6 ~ 20%. 

4. Relativistic corrections. This source of uncertainty appears because we treated r\b meson at the leading-order 
approximation in relative velocity. It can be estimated as ~ u 2 b ~ 10%. 

5. The uncertainties in the values of constants HI 6}) . For the three processes n^ — » J/i/jJ/t/jjJ/ipijj'ji^'tf)' these 
errors are estimated to be ~ 16%, 27%, 49%, respectively. 

6. Higher Fock states. It can be argued that at the scale [i relevant to r\b decay process, only a small fraction of 
quarkonium momentum is carried by the quark-gluon sea, typically ~ 5 — 10% [29(. Hence, we expect the effects 
of higher Fock states to be negligible, compared to other uncertainties considered here. 

The overall uncertainties of our calculations were obtained by adding the above errors in quadrature. 



C. Results and discussion 

By substituting the expressions for DAs and the necessary constants into eqs. ^ and ©, we get the following 
values for the three branching ratios: 

Br(j] h -> J/?pJ/ip) = (6.2 ± 3.5) x 10~ 7 , 
Br(r) b -> J/W) = (10±6)xl0~ 7 , (18) 
Br(rib-^ip'i>') = (3.7 ±2.8) x 10~ 7 . 
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It is interesting to compare these results with previous calculations. In particular, within the leading order NRQCD, 
one has ||: 

Br( Vb -» J/rl>J/rl>) = (2.4j£|) x 10~ 8 . (19) 

which is roughly 20 times smaller than our result shown above. The reason of this suppression can be traced to the 
expression for the amplitude ([9|), where all terms are in fact proportional to the constants f a and fs, which, in turn, 
are expressed through fv and fr (see eq. (JSJ). In the absence of relativistic and radiative corrections, the fine-tuning 
between fv, /t and the masses, clearly visible in eqs. (|12[) . guerantees that /a, fs and hence the formfactor F are 
proportional to (v 2 ), which is small for nonrelativistic systems. Taking relativistic and leading logarithmic radiative 
corrections to the constants f a and fs into account breaks the fine tuning, thus leading to a considerable enhancement 
of the branching ratio. To illustrate the above argument numerically, we take an infinitely narrow approximation for 
the DAs, parameters with fine-tuning given by eqs. fT2"]) , and (v 2 ) — 0.25, to obtain Br(i]b — > J/ipJ/tp) ~ 2 x 10 -8 , 
in agreement the leading order NRQCD result Q. Next, we take into account relativisitic and leading logarithmic 
radiative corrections to the constants j a and fs, but still use an infinitely narrow approximation for the DAs. In this 
case fine-tuning is broken, and we get ~ 3 x 10~ 7 , and order-of-magnitude increase compared to the NRQCD value. 
By including renormalization group evolution and relativistic motion into the DAs, we get a further increase of the 
branching ratio by a factor ~ 2. 

In Q the authors took into account one-loop radiative corrections and obtained 

Br(r) b -> J/ipJ/ip) = (2.1 - 18.6) x 10~ 8 . (20) 



Although this number seems to be compatible with ours shown in eq. (|18|) . we do not believe that the two results 
are in agreement with each other. In particular, the analytical form of the formfactor F obtained in [(| contains 
logarithmic terms: 



ReF ~ — log 



^log 2 ^ 
32 g M 2 



v/V 
1 Q M 2 

ImF ^16 l0g M^ + - < 21 > 

In the LC approach used in our calculation, all double logarithms cancel as the final partciles are colourless objects [30j . 
Moreover, there are only two reasons why a general QCD amplitude may contain large logarithms: renormalization 
and collinear divergences [30l [3l|. Clearly, the imaginary part of F is not renormalized at one loop, hence the large 
logarithm in eq. (|2ip must be due to a collinear divergence. However, it is known that collinear divergences can be 
factored out, and do not have an imaginary part [3l| . In light of these arguments, the result obtained in [|| looks 
strange. 

The authors of Q believe that there is no need for renormalization in their calculation of the radiative corrections, 
since the counterterms are proportional to the leading order contribution, which vanishes at the leading order in both 
a s and v c . We do not think that this statement is correct, since the expansion is done in operators which are not 
multiplicatively renormalizable. Therefore, the ultraviolet divergences may arise at the leading order in v c due to the 
i> c -suppressed operators. This effect violates NRQCD velocity scaling rules, and is discussed in detail in [20l l24j. 

Yet another estimate for the same branching ratio was obtained in Q , where the final-state interaction effects due 
to a different decay mechanism were taken into account, yielding 

Br(r] b -> J/ipJ/ip) = (0.5 x 10~ 8 - 1.2 x 10~ 5 ). (22) 

In conclusion, we have calculated the branching fractions of the decays r) b — > J ftp J /if), J/ipip', tp'tp' in the framework 
of the light cone formalism. The uncertainties of our calculation have also been assessed. Our results, presented in 
eqs. (p~8|) . are more than an order of magnitude larger than those obtained within NRQCD. 



Acknowledgments 

The authors thank A.K. Likhoded and A.V. Luchinsky for useful discussion. This work was partially supported by 
Russian Foundation of Basic Research under grant 07-02-00417. 



[1] B. Aubert et al. [BABAR Collaboration], Phys. Rev. Lett. 101, 071801 (2008) [Erratum-ibid. 102, 029901 (2009)] 
|arXiv:0807.1u86l [hep-ex]]. 



7 



E. Braaten, S. Fleming and A. K. Leibovich, Phys. Rev. D 63, 094006 (2001) [arXiv:hep-ph /0008091 . 
V. G. Kartvelishvili and A. K. Likhoded, Yad. Fiz. 40, 1273 (1984). 

G. T. Bodwin, E. Bra aten and G. P. Lepage, Phys. Rev. D 51, 1125 (1995) [Erratum-ibid. D 55, 5853 (1997)] 
|arXiv:hep-ph/9407339] . 

Y. Jia, Phys. Rev. D 78, 054003 (2008) [arXiv:hep-ph/0611130] . 

B. Gong, Y. Jia and J. X. Wang, Phys. Lett. B 670, 350 (2009) |arXiv:0808. 10341 [hep-ph]]. 
P. Santorelli, Phys. Rev. D 77, 074012 (2008) [arXiv:hep-ph/0703232l . 

E. Braaten and J. Lee, Phys. Rev. D 67, 054007 (2003) [arXiv:hep-ph/0211085]; 

K. Y. Liu, Z. G. He and K. T. Chao, Phys. Lett. B 557, 45 (2003) [arXiv:hep-ph/021118l]; 
K. Y. Liu, Z. G. He and K. T. Chao, Phys. Rev. D 77, 014002 (2008) |arXiv:hep-ph/0408141| . 



Y. J. Zhang, Y. j. Gao and K. T. Chao, Phys. Rev. Lett. 96, 092001 (2006) arXiv:hep-ph/0506076 
B. Gong and J. X. Wang, Phys. Rev. D 77, 054028 (2008) [arXiv:0712.4220l [hep-phjj. 
Y. J. Zhang, Y. Q. Ma and K. T. Chao, Phys. Rev. D 78, 054006 (2008) arXiv:0802.3655 [hep-ph]]. 
A. E. Bondar and V. L. Chernyak, Phys. Lett. B 612, 215 (2005) [arXiv:hep-ph/0412335[ 



V. V. Braguta, A K. Likhoded and A. V . Luchinsky, Phys. Rev. D 72, 074019 (2005) arXiv:hep-ph/0507275 

A. V. Berezhnoy, |arXiv:hep-ph/0703143| 

D. Ebert, R. N. Faustov, V. O. Galkin and A. P. Martvnenko. larXiv:0803.2124l [hep-ph], 
Z. G. He, Y. Fan and K. T. Chao, Phys. Rev. D 75, 074011 (2007) arXiv:hep-ph/0702239 . 
G. T. Bodwin, J. Lee and C. Yu , Phys. Rev. D 77, 094018 (2008) [arXiv:0710 .0995 [hep-ph]]. 
V. V. Braguta. larXiv:0811. 26401 [hep-ph], 

V. L. Chernyak and A. R. Zhitnitsky, Phys. Rept. 11 2, 173 (1984). 

V. V. Braguta, A. K. Likhoded and A. V. Luchinsky. [arXiv:0902.0459l [hep-ph]. 



P. Ball, V. M. Braun, Y. Koike and K. Tanaka, Nucl. Phys. B 529, 323 (1998) arXiv:hep-ph/9802299 
V. V. Braguta, A. K. Likhoded and A. V. Luchi nsky, Phys. Lett. B 646, 80 (2007) |arXiv7h ep-ph/0611021 
V. V. Braguta, Phys. Rev. D 75, 094016 (2007) [arXiv:hep-ph/0701234| . 
V. V. Braguta, Phys. Rev. D 77, 034026 (2008) arXiv:0709. 3885 hep-ph]]. 
V. V. Braguta, A. K. Likhoded and A. V. Luchinsky, arXiv:0810.3607 [hep-ph]. 
V. V. Braguta, Phys. Rev. D 78, 054025 (2008) [arXiv:0712. 14751 [hep-ph]] . 
V. G. Kartvelishvili and A. K. Likhoded, Sov. J. Nucl. Phys. 42 (1985) 823. 
G. P. Lepage and S. J. Brodsky, Phys. Rev. D 22, 2157 (1980). 
A. V. Smilga and M. I. Vysotsky, Nucl. Phys. B 150, 173 (1979). 



